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NOTATION 


Airfoil  coefficients  (used  as  design  variables) 

Chordwise  force  coefficient 

Drag  coefficient  (total) 

Lift  coefficient 

Moment  coefficient 

Normal  force  coefficient 

Pressure  coefficient 

Momentum  coefficient 

Chord  length 
Objective  function 
Constraint  function 


Slot  height 

Length  of  a  panel 

Mach  number 

Number  of  panels 

Source  strength 

Direction  of  search 

Arc  length  along  airfoil  surface 

Thickness 

Velocity 


Vector  consists  of  design  variables 


Airfoil  coordinates 
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a  Angle  of  attack 

* 

a  Scalar  defining  distance  of  travel 

T  Total  circulation 

Y  Vortex  strength 

* 

6  Boundary  layer  displacement  thickness 

p  Density 

Xj  Design  variables 

Superscripts 
l  Lower 

m  Iteration  number  during  optimization 

u  Upper 

Subscripts 

e  Edge  of  boundary  layer 

f  Frictional 

i  ith  panel 

j  jth  design  variable;  also  jet 

m  Iteration  number  during  analysis 

max  Maximum 

min  Minimum 

p  Pressure 

sepl  Separation  point,  lower  surface 

sepu  Separation  point,  upper  surface 
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ABSTRACT 


A  numerical  procedure  developed  for  optimizing  the 
circulation  control  airfoils  is  presented.  The  pro¬ 
cedure  finds  the  optimum  basic  airfoil  shapes  subjected 
to  specified  flow  conditions  and  geometric  constraints. 

It  consists  of  a  numerical  optimization  code  for  linear 
or  nonlinear  constrained  problems  coupled  with  a  viscous- 
potential  flow  interaction  analysis  for  necessary 
viscous-inviscid  flow  field  calculations.  The  desired 
airfoil  shape  is  defined  by  a  combination  of  baseline 
shapes  representative  of  airfoils  suitable  for  circu¬ 
lation  control  purposes.  The  coefficients  of  these 
basis  vectors  are  then  used  as  the  design  variables  in 
the  optimization  process.  Three  baseline  shapes  (a 
cambered  ellipse,  a  cambered  ellipse  with  a  drooped 
trailing  edge,  and  a  cambered  ellipse  with  a 
logarithmically  spiralled  trailing  edge)  are  employed 
for  special  contouring  of  the  trailing  edge  geometry. 

With  some  minor  modification  of  the  analysis  method, 
the  combined  program  allows  optimization  for  maximum 
lift  without  substantial  difficulty;  but  for  minimizing 
the  drag,  further  improvement  of  the  analysis  method  is 
required. 

Numerical  results  were  obtained  for  maximizing  CL 

with  a  blowing  coefficient  of  0.04,  a  Mach  number  of 

0.3,  and  at  angles  of  attack  of  -5  deg  and  -2  deg. 
Contraints  on  the  lift  augmentation  ratio  and  the  lift-to- 
drag  ratio  were  imposed.  The  lift  coefficient  increased 
from  an  initial  value  of  2.05  to  a  final  value  of  2. 34 
for  the  case  of  a  =  -5  deg.  A  similar  result  was  ob¬ 
tained  for  the  case  of  a  =  -2  deg  with  the  same  flow 
condition.  The  resulting  airfoil  profile  is  somewhat 
between  the  baseline  shapes  of  cambered  ellipse  and  the 
drooped  trailing  edge.  The  result  agrees  qualitatively 
with  the  available  experimental  data  for  which  the  air¬ 
foils  having  cambered  elliptic  and  drooped  trailing  edge 
contours  exhibit  better  performance  than  the  one  with 
spiralled  trailing  edge  geometry.  Efforts  for  improving 
the  drag  prediction  capability  of  the  analysis  method 
are  in  progress. 
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INTRODUCTION 

In  the  development  of  technology  for  V/STOL  aircraft,  efforts  have 
been  made  to  increase  the  lift  coefficient  through  blowing.  Systems 
which  use  blowing  may  be  classified  into  three  major  categories: 

(1)  the  blowing  flap,  (2)  the  jet  flap,  and  (3)  the  circulation  control 
airfoils.  In  the  application  of  blowing  to  helicopter  rotors,  the 
circulation  control  (CC)  system  offers  two  basic  advantages,  namely,  the 
much  higher  lift-to-thrust  augmentation  for  a  given  slot  momentum  and  the 
mechanical  simplifications  in  lift  control. 

Experimental  data  are  available  for  circulation  control  air- 
2-4* ** 

foils.  ’  Analytical  methods  developed  earlier  used  an  integral 

approach  to  model  the  flow  as  an  incompressible  turbulent  boundary  layer 
flow  mixed  with  a  wall  jet.^’^  Later,  Dvorak  and  Kind  solved  the  wall 
jet  flow  region  problem  by  the  finite  difference  scheme  but  retained  the 
integral  approach  for  flow  regions  where  ordinary  boundary  layers  prevail.^ 
Reasonably  good  lift-momentum  results  have  been  reported  in  using  these 
methods,  although  a  major  deficiency  in  drag  prediction  still  exists. 
Attempts  to  incorporate  a  more  realistic  eddy  viscosity  model  are  found 
elsewhere.^ 

To  meet  the  performance  requirement  for  several  naval  applications,*** 

an  optimization  procedure  is  highly  desirable  to  provide  an  advanced 

design  capability.  In  the  present  work,  a  numerical  approach  is  taken  to 

optimize  the  airfoil  for  maximum  lift  or  for  minimum  drag,  or  for  maximum 

lift-to-drag  ratio.  Various  constraints  including  thickness,  lift-to-drag 

ratio,  and  lift-momentum  flux  relations  are  imposed  to  generate  optimal 

airfoil  shapes.  In  so  doing,  a  numerical  optimization  scheme  of 
9 

Vanderplaats  for  linear  or  nonlinear  constrained  problems  has  been  em¬ 
ployed  and  coupled  with  a  viscous-potential  flow  interaction  analysis  of 
Dvorak  and  Kind^  for  necessary  viscous-inviscid  flow  field  calculations. 

*A  complete  listing  of  references  is  given  on  page  39. 

**Work  performed  by  J.S.  Abramson  for  two-dimensional  subsonic  wind 
tunnel  tests  of  the  circulation  control  airfoils  documented  as  DTNSRDC  TM- 
16-76-42  of  December  1975. 

***During  1977  to  1979,  the  circulation  control  airfoils  have  been 
applied  to  an  A-6  testbed  aircraft,  an  H-2  helicopter,  and  a  stopped  rotor 
(X-Wing)  aircraft  at  DTNSRDC. 
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Reasonably  good  agreement  between  theoretical  and  experimental  results 
for  certain  cases  has  been  reported  to  warrant  some  confidence  in  the 
method . 

ANALYSIS  METHOD 

A  schematic  of  the  flow  about  a  CC  airfoil  is  shown  in  Figure  1.  As 
opposed  to  the  conventional  airfoil,  a  typical  CC  airfoil  is  equipped 
with  a  blowing  slot  on  the  upper  surface  for  energizing  the  flow  in  the 
viscous  layer  and  a  rounded  trailing  edge  for  deflecting  the  jet.  High 
lift  can  be  generated  because  of  increased  circulation  created  by  the 
blowing  jet.  The  flow  is  characterized  by  outer  inviscid  flow  and  inner 
viscous  flow  consisting  of  boundary  layers,  wall  jet,  and  separation 
bubble. 

Rather  than  starting  from  scratch,  the  analysis  method  developed 
by  Dvorak  and  Kind^  with  some  minor  modifications  was  adopted.  The 
method  is  found  to  be  the  most  comprehensive  one  available.  A  brief 
description  of  the  theory  is  given  below. 

DVORAK-KIND  METHOD 

Dvorak  and  Kind^  considered  the  subject  problem  as  a  viscous-inviscid 
interaction  flow  over  a  two-dimensional  airfoil  section  with  or  without 
slot  blowing  on  the  upper  surface.  The  inviscid  potential  flow  is  calcu¬ 
lated  first  and  then  followed  by  boundary  layer  and  wall  jet  developments 
which  are  computed  using  pressure  distributions  obtained  from  the  potential 
flow  analysis.  The  viscous  effects  are  modeled  and  the  procedure  is 
repeated  until  a  converged  solution  is  achieved. 

Potential  Flow  Solution 

The  potential  flow  is  represented  by  vorticity  distribution  along  the 
panels.  The  vortex  strength  is  assumed  to  vary  linearly  along  each  panel 
and  is  continuous  at  all  junction  points.  The  horizontal  and  vertical 
component  of  the  velocity  induced  by  the  vorticity  according  to  the 
potential  flow  theory  is 
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where  y  is  the  vortex  strength  along  the  panel.  The  solution  procedure  is 
similar  to  those  for  conventional  airfoils  except  that  the  Kutta  condition 
is  applied  differently.  For  a  rounded  trailing  edge,  which  is  typical  in 
a  CC  airfoil,  the  Kutta  condition  is  replaced  by  an  equation  which 
specifies  the  value  of  the  total  circulation  T  around  the  airfoil 
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where  2,^  is  the  length  of  the  ith  panel.  The  value  of  T,  together  with  the 
blowing  coefficient  C^,  has  direct  bearing  on  the  lift  coefficient.  A 
detailed  description  of  the  potential  flow  solution  can  be  found  in 
Reference  10. 


Laminar  and  Turbulent  Boundary  Layers 

The  forward  and  aft  stagnation  points  obtained  by  the  potential  flow 
solution  divide  the  flows  between  the  upper  and  lower  surfaces  of  the  air¬ 
foil.  Based  on  pressures  obtained  from  the  potential  flow  solution, 

calculation  of  boundary  layers  starts  with  the  Hiemenz  stagnation  flow 
11  12 

solution.  Curie's  integral  method  is  employed  for  calculating  the 

laminar  boundary  layers  downstream  of  the  stagnation  point  for  both  upper 

and  lower  surfaces.  After  the  transition  predicted  by  Granville's 
13 

empirical  formula,  the  turbulent  boundary  layers  for  the  remaining  lower 
surface  to  the  separation  point,  and  the  upper  surface  to  the  blowing  slot 
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are  calculated  using  the  Nash-Hicks  integral  method.  At  the  slot,  the 
turbulent  boundary  layer  mixes  with  a  wall  jet  of  a  uniform  velocity 
distribution.  A  finite  difference  scheme  is  used  in  solving  the  mixed  flow. 
The  flow  proceeds  around  a  highly  curved  surface  with  very  strong  adverse 
pressure  gradients.  The  effect  of  surface  curvature  downstream  of  the 
blowing  slot  is  accounted  for  by  including  the  normal  momentum  equation  for 
radial  variation  of  static  pressure  in  the  finite  difference  procedure. 

The  flow  eventually  separates  after  passing  the  blunted  trailing  edge  and  a 
separation  pressure  is  noted. 

Viscous-Potential  Flow  Interaction 

The  viscous  effects  produced  by  the  boundary  layer  development  are 
modeled  by  a  source  distribution  along  the  airfoil  surface.  The  source 
strength  q^  at  any  junction  point  of  a  surface  panel  is  obtained  by 
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where  s  =  arc  length 

Vg  =  local  external  velocity 
* 

5  =  boundary-layer  displacement  thickness 

Equation  (A)  applies  to  the  conventional  boundary  layers  as  well  as  to  the 

* 

wall  jet.  In  the  latter  case,  the  product  tends  to  decrease,  yielding 

a  nega  ive  source  strength  (i.e.,  sink). 

The  contribution  of  q^  is  then  implemented  in  the  potential  flow 
solution  which  alters  the  results  of  the  pressure  distribution.  The  entire 
procedure  is  then  repeated  with  subsequent  values  of  the  total  circulation 
estimated  by 
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where  the  numerical  constant  k  has  a  value  from  0.1  to  0.3.  The  calcu¬ 
lations  are  terminated  when  the  lift  coefficient,  which  is  directly 
related  to  the  value,  remains  within  a  specified  tolerance  and,  at  the 

same  time,  the  differences  in  pressure  coefficients  at  the  upper  and  lower 

ends  of  the  separation  bubble,  /C  -C  \  diminish. 

\  psepu  Psepl  / 

The  lift  and  pressure  drag  coefficients  are  then  evaluated  by  inte¬ 
grating  the  surface  pressure  distribution  using  the  trapezoidal  rule 
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C.,  =  C  cos  a  -  C  sin  a;  C_.  =  C  sin  a  +  C  cos  a 

L  n  c  D  n  c 

F 

where  x  and  z  are  airfoil  coordinates  normalized  by  the  chord  length,  and 
N  is  the  number  of  panels.  The  total  drag  is  the  sum  of  the  pressure  drag 
and  skin  friction,  less  the  blowing  momentum  coefficient 
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(7) 


The  momentum  coefficient  is  defined  by 


:  =  ,Pj  Vj  h 
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where  PjVj  is  the  momentum  flux  of  the  blowing  jet  and  h  is  the  jet  height 
at  the  slot. 
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COMPUTER  PROGRAM 

The  method  has  been  coded  In  FORTRAN  and  is  known  as  the  CIRCON 
program.  It  has  been  written  in  an  overlay  form  to  reduce  core  require¬ 
ments.  A  detailed  description  of  the  code  can  be  found  in  Reference  15. 

COMPARISON  BETWEEN  THEORY  AND  EXPERIMENT 

Previous  results  of  the  Dvorak-Kind  method,  reported  in  Reference  7 
show  good  agreement  between  the  theory  and  the  experiment.  Cases  con¬ 
sidered  in  Reference  7  include  comparison  of  the  computed  results  with 
experimental  data  of  a  cambered  20-percent  ellipse. 

The  method  has  been  applied  to  other  cases  at  DTNSRDC  on  different 
airfoils  with  various  flow  conditions,  among  which  is  discussed  the  case 
of  a  15-percent  cambered  ellipse  at  M^  =  0.3  and  a  =  -0.01  deg.  The 
computed  results  are  compared  with  experimental  data*  obtained  at  the 
DTNSRDC  7  ft  x  10  ft  transonic  wind  tunnel  where  the  blockage  effects  can 
be  considered  insignificant  at  moderate  blowing  conditions.  The  results 
are  given  in  Table  1.  Figures  2  and  3  show  the  calculated  lift  and  drag 
coefficients  at  various  blowing  conditions  between  C^  =  0.0073  and 
0.0384,  along  with  experimental  data.  Qualitative  agreement  between  the 
theoretical  and  experimental  lift  coefficients  is  obtained.  However, 
serious  discrepancies  in  drag  coefficients  are  observed.  In  particular, 
the  dip  in  the  theoretical  drag  curve  is  not  detected  in  the  experiment. 
Some  representative  pressure  distributions  (for  C^  =  0.0111,  0.0160,  and 
0.0220)  are  given  in  Figure  4,  which  may  reveal  the  disagreement. 

The  results  indicate  that  although  the  Dvorak-Kind  method  is  the 
most  sophisticated  approach  available  for  analyzing  the  flow  over  a  CC 
airfoil,  it  is  still  inadequate.  Efforts  in  numerical  optimization  here, 
discussed  later,  will  be  restricted  to  the  cases  using  the  lift  co¬ 
efficient  as  the  objective  function  to  obtain  optimal  shapes.  Nonetheless, 
the  approach  is  capable  of  achieving  qualitative  results  which  may  be  used 
to  guide  the  design  engineers. 


*Experimental  data  were  provided  by  J.B.  Wilkerson  and  J.S.  Abramson 
of  DTNSRDC. 
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OPTIMIZATION  PROGRAM 

The  aforementioned  analysis  method  has  been  coupled  with  a  numerical 

9 

optimization  scheme  developed  by  Vanderplaats.  A  brief  description  of 
the  optimization  technique  is  given  here  and  the  procedure  for  automating 
the  CC  airfoil  design  is  outlined. 

VANDERPLAATS  METHOD 

The  optimization  program  is  based  on  the  method  of  feasible  directions 
for  constrained  problems.  The  optimization  problem  is  stated  mathematical¬ 
ly  as 


Minimize  OBJ  =  F(X) 


which  is  subject  to: 


G .  (X)  <  0 
J  “ 


(j=l,m) 


(9) 


l  .  .  u 

X,  1  X  1  X 


(10) 


where  OBJ  is  the  objective  function.  The  vector  X  contains  n  design 
variables.  The  function  G^ (X)  defines  the  constraints  which  the  designer 
wishes  to  impose  on  the  optimization  problem.  Functions  F(X)  and  G^ (X) 
may  be  either  implicit  or  explicit  of  the  design  variables  X  but  must  be 

&  VI 

continuous.  Variables  Xj  and  Xj  define  the  lower  and  upper  bounds, 
respectively,  on  the  design  variable  and  are  the  limits  over  which  F(X) 
and  Gj (X)  are  defined.  If  the  inequality  condition  of  Equation  (9)  is 
violated  (G^ (X)  >  0)  for  any  constraint,  the  constraint  is  said  to  be 
violated.  If  the  equality  condition  is  met  (G^ (X)  =  0),  the  constraint  is 
called  active,  and  if  the  strict  inequality  is  met  (Gj (X)  <  0),  the 
constraint  is  inactive.  Because  a  precise  zero  is  seldom  meaningful  in 
the  digital  computation,  a  constraint  is  called  active  if  its  value  is 
within  a  specified  tolerance. 
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The  n-dimensional  space  spanned  by  the  design  variables  X  is  referred 
to  as  the  design  space.  Any  design  which  satisfies  the  inequalities  of 
Equations  (9)  and  (10)  is  referred  to  as  a  feasible  design.  If  a  design 
violates  ont  jr  more  of  these  inequalities,  it  is  said  to  be  infeasible. 

The  minimum  feasible  design  is  said  to  be  optimal.  Note  that  if  one  wishes 
to  maximize  some  function  such  as  lift,  it  can  be  done  simply  by  minimizing 
the  negative  of  lift.  Thus,  any  design  problem  can  be  cast  in  the  fore¬ 
going  form. 

The  optimization  program  begins  with  an  initial  X  vector  which  is 
input  to  the  program  and  may  or  may  not  define  a  feasible  design.  The 
optimization  process  then  proceeds  iteratively  by  the  following  recursive 
relationship 


_ m+1  _  .  m 

X  =  X  +  a  S  (11) 

where  m  =  iteration  number 

S  (vector)  =  direction  of  search  in  the  n-dimensional 
design  space 

* 

a  =  scalar  which  defines  the  distance  of 

travel  in  direction  S 
_m 

The  S  is  obtained  by  moving  in  the  direction  of  the  steepest  descent  (the 

negative  gradient  of  the  objective  function)  without  violating  constraints 
* 

The  scalar  a  is  determined  by  a  one-variable  search  based  on  a  polynomial 
fit  of  several  trial  values. 

The  procedure  is  repeated  with  the  aid  of  a  conjugate  direction 
algorithm^  in  determining  the  new  search  direction.  When  the  constraint 
is  encountered  in  the  searching  process,  the  new  search  direction  is  found 
using  Zoutendijk's  method  of  feasible  directions. ^  The  optimum  point  is 
achieved  where  no  direction  can  be  found  that  will  reduce  the  objective 
without  violating  the  constraints. 

COMPUTER  PROGRAM 

q 

The  optimization  procedure  is  coded  in  FORTRAN  known  as  CONMIN.  The 
program  has  been  expanded  since  its  first  appearance  in  1973.  A  main 
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control  program  known  as  COPES  having  four  modes  of  application,  namely, 
ANALYSIS,  SENSITIVITY  STUDY,  TWO-VARIABLE  SPACE,  and  OPTIMIZATION  is  now 
available.  Also  an  approximation  optimization  procedure  is  added  using  a 
truncated  Taylor’s  series  to  reduce  the  number  of  required  analyses. 

Both  old  and  new  versions  have  been  widely  used  in  various  engineering 
disciplines. 

OPTIMIZATION  OF  CIRCULATION  CONTROL  AIRFOIL 
The  subject  problem  is  to  optimize  the  airfoil  for  maximum  lift,  or 
for  minimum  drag,  or  for  maximum  lift-to-drag  ratio.  Various  constraints 
including  minimum  lift  range,  maximum  drag  range,  angle  of  attack,  lift 
augmentation  ratio,  thickness  ratio,  jet  detachment,  and  linear  lift- 
momentum  flux  relations  are  to  be  imposed  to  generate  optimal  airfoil 
shapes. 

The  optimization  model  for  the  problem  contains  the  objective,  the 
constraints,  and  the  design  variables.  For  example,  if  we  want  to 
optimize  the  airfoil  for  maximum  lift  subject  to  thickness  ratio  and  lift- 
to-drag  ratio  constraints,  the  optimization  model  would  be 

Objective  -  F(X):  -  CL 

Constraints  -  G(X):  (t/c)  <  t/c  <  (t/c)  (12) 

min  —  —  max 

<Cl/S>  ±  CL/CD  i  <VCD>  (13) 

min  max 

Design  Variables  -  X:  a^ ,  j  =  1,  n 

where  the  variables  a^  are  the  coefficients  of  a  linear  combination  of  n 
baseline  airfoils  that  represent  the  design  airfoil: 


The  lift  coefficient  and  the  drag  coefficient  are  nonlinear  implicit 
functions  of  the  design  variables  to  be  obtained  by  the  analysis  method 
and  the  thickness  ratio  is  a  linear  function  of  the  design  variables 
determined  directly  by  Equation  (14). 

The  model  is  implemented  numerically  by  interaction  between  the 
optimization  and  airfoil  analysis  programs.  The  program  flow  chart  for 
such  an  operation  is  shown  in  Figure  5.  In  so  doing,  the  airfoil  analysis 
program  CIRCON  is  merged  with  the  optimization  code  CONMIN  on  the  CDC  7600 
computer  at  the  NASA  Ames  Research  Center  to  form  an  airfoil  design 
capability.  The  overlay  of  the  combined  program  is  shown  in  Figure  6. 

The  program  can  be  operated  either  at  Ames  or  at  DTNSRDC  using  a  remote 
terminal. 

AIRFOIL  REPRESENTATION 
Trailing  Edge  Contouring 

The  trailing  edge  geometry  has  direct  bearing  on  the  jet  deflection 
characteristics  and,  therefore,  our  first  effort  is  to  contour  the  blunt 
trailing  edge  of  the  CC  airfoil.  Three  baseline  shapes,  i.e.,  those  of  a 
cambered  ellipse,  a  cambered  ellipse  with  a  drooped  trailing  edge,  and  a 
cambered  ellipse  with  a  logarithmically  spiralled  trailing  edge,*  are 
employed  in  test  cases.  These  shapes  (shown  in  Figure  7),  representative 
of  airfoils  suitable  for  circulation  control  purposes,  allow  special 
contouring  of  the  trailing  edge  geometry  which  directly  affects  the  air¬ 
foil  performance.  Equation  (10),  therefore,  simplifies  to 

Y  =  ax  Yx  +  a2  Y2  +  a3  Y3  (15) 

The  exact  coordinates  of  these  baseline  airfoils  are  listed  in  Table  2. 

All  three  shapes  have  identical  coordinates  up  to  x/c  =  0.95. 
Variation  is  allowed  between  0.95  x/c  1.00.  For  this  special  case, 
since  the  thickness  of  the  resulting  airfoil  remains  constant,  we  have 

*These  baseline  shapes  are  provided  by  E.O.  Rogers  of  the  Aviation 
and  Surface  Effects  Department,  DTNSRDC. 
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al  +  a2  +  a3  =  1 


(16) 


Accordingly,  Equation  (15)  can  be  written  as 

Y  =  ax  Yx  +  a2  Y2  +  (l-a^)  Y3  (17) 

The  new  form  of  Equation  (17)  reduces  the  actual  number  of  design 
variables  from  three  to  two  (but  still  three  baseline  shapes).  It,  in 
turn,  reduces  the  number  of  gradients  and,  thus,  simplifies  the  search 
process.  It  also  eliminates  the  thickness  constraint,  which  is  auto¬ 
matically  satisfied.  Note  that  since  negative  values  of  a^,  a^,  or  a^  are 
allowed,  the  resulting  airfoil  shape  can  be  drastically  different  from 
these  original  profiles. 

Overall  Airfoil  Shaping 

The  ultimate  purpose  of  the  present  work  is  to  derive  a  CC  airfoil 
for  optimum  performance.  The  work  involves  overall  airfoil  shaping  based 
on  more  general  representation  of  baseline  profiles.  The  computer  program 
has  been  set  up  to  accept  as  many  as  twenty-five  (25)  baseline  shapes, 
although  practically  six  would  be  sufficient  to  cover  a  broad  range  of 
interest. 

MODIFICATION  OF  ANALYSIS  PROGRAM 

A  major  difficulty  in  using  the  viscous- invisc id  interaction  method 
together  with  the  optimization  code,  however,  is  its  inability  in 
providing  fairly  smooth  gradients  of  the  objective  function  when  the 
design  variables  are  perturbed  in  the  course  of  optimization.  The  problem 
is  a  common  one  in  most  viscous-inviscid  interaction  methods  in  which  the 
viscous  and  inviscid  flows  are  calculated  separately  and  the  final  solution 
is  reached  by  the  iterative  process.  In  addition,  as  mentioned  earlier, 
the  correlation  between  the  analysis  method  and  the  experiment  requires 
further  improvement.  The  problem  is,  therefore,  twofold,  i.e.,  con¬ 
vergence  and  accuracy. 
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To  circumvent  the  problem,  some  modifications  to  the  analysis  program 
need  to  be  carried  out.  The  approach  is  to  make  the  program  operational 
with  the  least  modification  required  and  then  to  refine  it  as  necessary. 

The  first  step  is,  therefore,  to  smooth  out  some  irregularities  in 
aerodynamic  force  coefficients,  which  might  be  partially  caused  by  the 
simple  trapezoidal  rule  used  in  integrating  the  surface  pressure  distri¬ 
bution  on  a  limited  number  of  panels.  Equation  (6)  is  thus  replaced  by  a 
second-order  Lagrangian  interpolation  scheme  for  the  integrand,  which 
results  in 


N+l  _ 

„  y  c  xi+i  xj-i 
■n  Z*  Pi  x  -x 
i=l  L  1 


+2  +  C 


xj-xj-l 
l+l  l-l 


(xi+i~xi)2  1  (xi+rxi)  nR. 

L-l  _ (xi~xi-l) (xi+l-xi-l) .  6 


For  evaluating  Cc,  x's  are  replaced  by  z's. 

The  difference  between  the  two  schemes  can  be  appreciable  in  drag 
values  attributed  to  the  leading  and  trailing  edge  regions  where  change  in 
sign  of  the  pressure  coefficient  takes  place.  Figures  8a  and  8b  indicate 
the  values  (the  net  areas,  which  are  the  differences  of  the  shaded 
regions)  resulting  from  these  two  schemes.  The  advantage  of  using  the 
Lagrangian  scheme  is  twofold:  (a)  improving  the  convergences  of  values 
and  (b)  improving  the  accuracy  for  both  CL  and  C^.  This  is  shown  in 
Figure  9.  The  slope  of  the  CT  variation  during  the  iteration  for  the 
Lagrangian  scheme  is  smaller  than  that  for  the  trapezoidal  rule.  After 
eight  iterations  by  using  Equation  (18),  the  resulting  (converged  since 
AC  /C  <  e  where  e  is  a  specified  tolerance)  is  smaller,  but  C  is  larger, 
than  it  would  be  if  using  Equation  (6).  The  accuracy  of  both  C  and  C  has 
been  improved  since  the  theory  overpredicts  the  lift  but  underpredicts  the 
drag  (see  Figures  2  and  3).  It  is  anticipated  that  if  a  large  number  of 
panels  is  used,  the  difference  between  the  two  schemes  should  be  minimal. 
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The  modification  (although  a  minor  one),  when  coupled  with  the  recast 
equation  for  airfoil  representation  (Equation  (17)),  allows  performing 
optimization  runs  for  maximizing  the  lift  without  substantial  difficulty; 
but  for  minimizing  the  drag,  further  improvement  of  the  analysis  method  is 
required.  This  will  be  discussed  further  in  the  next  section. 

The  second  phase  of  the  effort,  which  is  now  in  progress,  will  be 
concentrated  on  modifications  of  the  analysis  program  including  (a)  remodel¬ 
ing  the  potential  flow  solution  to  consider  the  effect  of  the  large  wake 
associated  with  the  blowing  jet  and  (b)  removal  of  the  assumption  of 
constant  pressure  inside  the  separation  bubble  which  has  caused  some 
arbitrary  variation  in  aerodynamic  force  coefficients.  Although  only 
slightly  noticeable  in  lift  coefficients,  these  variations  produce  wide¬ 
spread  differences  in  drag  values.  Finally,  (c)  rearrangement  of  the 
viscous-potential  flow  interaction  procedure  to  improve  convergence  of  the 
overall  solution. 


RESULTS  AMD  DISCUSSION 

Numerical  results  were  calculated  on  the  NASA  Ames  CDC  7600  computer 
using  a  remote  terminal  at  DTNSRDC.  The  merged  program  was  first  checked 
by  the  ANALYSIS,  SENSITIVITY  STUDY,  and  TWO-VARIABLE  SPACE  runs  in 
accordance  with  the  control  options  of  the  main  program  COPES  before  ex¬ 
tensive  OPTIMIZATION  runs  were  performed.  Results  of  the  ANALYSIS  runs  are 
tabulated  in  Table  1,  plotted  in  Figures  2,  3,  and  4.  and  have  been 
discussed  previously  in  the  Analysis  Method  section.  Results  of  other 
typical  runs  are  summarized  in  Table  3  and  presented  in  Figures  10  through 
16. 

The  first  optimization  run  was  a  simple  case  for  minimizing  the 
frictional  drag  with  the  angle  of  attack  as  the  design  variable.  The 
results  of  this  run  yield  a  slight  decrease  in  skin  friction  correspond¬ 
ing  to  a  change  of  the  design  variable  a  from  -5  deg  to  -4.9  deg.  The 
purpose  of  the  run  was  to  test  the  status  of  the  program  rather  than 
actually  design  for  minimum  skin  friction,  which  is  just  a  small  portion 
of  the  total  drag  and  has  only  a  small  range  of  variation. 
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The  next  two  runs  were  aimed  at  minimizing  the  total  drag  subject  to 
constraints  of  C  /C  ,  Cn  /C  ,  and  t/c  values  (Runs  6  and  7  in  Table  3). 

l*  U  L* 

The  angle  of  attack  was  used  as  the  design  variable  for  Run  6,  and  the 
coefficients  a^,  a and  a^  for  Run  7.  Unfortunately,  there  was  no 
change  in  the  objective  function  for  either  case. 

Since  the  skin  friction  values  remain  fairly  constant,  the  variation 
of  the  total  drag  is  basically  attributed  to  the  variation  of  the  pressure 
drag.  From  Runs  8  through  10,  therefore,  the  pressure  drag,  which  is  more 
direct  than  the  total  drag,  was  employed  as  the  objective  function  under 
similar  constraints.  With  the  aid  of  the  recasted  airfoil  equation  and 
the  modifications  in  the  integration  scheme  of  the  analysis  method, 
a  moderate  change  in  the  pressure  drag  is  obtained  in  Run  10b,  as  shown  in 
Figure  10.  The  advantage  of  using  the  Lagrangian  scheme  is  clearly 
demonstrated  in  Figure  10  although  the  general  outcome  in  these  cases 
appears  not  as  good  as  anticipated. 

The  problem  in  minimizing  the  drag  can  best  be  understood  by  the  plot 

of  the  CD  value  in  two-variable  space  given  in  Figure  11.  The  zigzag  CD 

P  P 

distribution  throughout  the  a^  -  a^  space  cannot  provide  smooth  values  of 

gradients  9CQ  /9a^,  and  90^  /da^  and,  therefore,  the  optimization  technique 
P  P 

fails  to  respond  properly  when  a^  and  a2  are  perturbed.  There  are,  how¬ 
ever,  still  some  regions  where  local  continuity  of  the  gradients  exists. 
These  permit  changes  in  the  objective  function  during  the  course  of 
optimization,  such  as  those  shown  in  Figure  10.  Nevertheless,  the  result 
of  Figure  10  can  only  be  regarded  at  best,  as  locally  optimal;  its  useful¬ 
ness  is  very  limited. 

Subsequent  efforts  were  concentrated  on  cases  using  the  lift  co¬ 
efficient  as  the  objective.  Numerical  results  were  obtained  for  maximizing 
with  a  blowing  coefficient  C  of  0.04,  at  a  Mach  number  of  0.3,  and  at 
angles  of  attack  of  -5  deg  and  -2  deg.  Constraints  of  the  lift  augmen¬ 
tation  ratio,  50  £  C  /C  _<  60,  and  the  lift-to-drag  ratio,  -0.014  <  C  /C 

Li  L)  Li 

_<  -0.006,  were  imposed.  Figure  12  shows  the  values  of  C  versus  the 

Li 

search  position  during  the  course  of  optimization  for  the  case  of 
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a  =  -5  deg.  The  search  starts  with  a^  =  0.3333  and  a ^  -  0.3334,  which 
yields  =  2.054.  The  search  stops  after  the  change  of  the  objective 
function  within  a  specified  tolerance,  consecutively  for  three  times.  A 
final  lift  coefficient  of  2.34,  is  achieved  (Run  12  in  Table  3)  with 
design  variables  a^  =  0.4084,  =  0.4798,  and  a^  =  0.1180.  A  similar  run 

with  less  restrictive  constraints,  but  using  three  design  variables  along 
with  the  original  trapezoidal  rule  integration  (Run  11a  in  Table  3),  is 
also  plotted  in  Figure  12.  Again,  it  indicates  that  the  advantage  of 
using  the  recast  airfoil  equations  along  with  the  Lagrangian  scheme  is 
distinctive. 

Figure  13  shows  the  C  values  in  the  two-variable  design  space. 

L 

Continuity  of  C  values  in  most  regions  is  maintained,  as  marked  by  solid 
lines.  Some  uncertainty  is  involved  in  the  region  with  broken  lines.  As 
a  double  check,  the  intermediate  C  values  of  Run  12  are  superimposed  in 

1j 

Figure  13.  Figure  13  indicates  that  the  optimal  value  approaches  the 
true  maximum  fairly  closely.  The  C  value  tends  toward  but  terminates 

L 

before  the  maximum  because  the  termination  criterion  has  been  met.  The 
test  warrants  further  extension  to  more  baseline  shapes  in  representing 
the  airfoil  so  that  greater  versatility  may  be  achieved. 

The  resulting  airfoil  profile,  which  is  shown  in  Figure  14,  is 
somewhat  between  baseline  shapes  and  The  result  agrees  quali¬ 

tatively  with  the  available  experimental  data  for  which  the  airfoils 
having  cambered  elliptic  and  drooped  trailing  edge  contours  exhibit  better 
performance  than  the  spiralled  one.  It  is  of  interest  to  note  that  at 
negative  angles  of  attack,  where  the  CC  airfoils  normally  operate,  high  CL 
values  are  usually  accompanied  by  low  drag  coefficients. 

Figure  15  shows  values  during  the  optimization  process  for  the 
case  of  =  0.3,  C  =  0.04,  and  a  =  -2  deg.  The  calculation  was 

OO  ° 

terminated  after  five  searches:  the  values  gained  from  2.255  to  2.740 
(Run  13  in  Table  3).  The  final  design  ended  up  with  a^  =  0.499, 

=  0.624,  and  a^  ■  -0.123.  The  resulting  airfoil  is  depicted  in 
Figure  16.  Note  that  a^  has  a  negative  value,  which  yields  a  more  drooped 
trailing  edge  than  the  previous  case  of  a  =  -5  deg. 


Finally,  the  effect  of  the  constraint  on  the  resulting  values 
is  examined.  This  is  shown  in  Figure  15,  where  the  results  of  Run  14 
(with  constraints  encountered  during  the  course  of  optimization)  are 
compared  with  those  of  Run  13  (without  constraints  encountered). 
Considerable  penalty  due  to  the  active  constraints  are  realized. 

CONCLUSIONS 

A  numerical  procedure  for  optimizing  circulation  control  airfoils  is 
developed.  Based  on  the  first  phase  of  work  in  the  continuing  effort  of 
improving  aerodynamic  performance  of  the  circulation  control  airfoil — the 
following  conclusions  may  be  drawn: 

1.  A  significant  gain  (about  12  to  15  percent)  in  lift  coefficients 
of  CC  airfoils  may  be  achieved  by  optimizing  the  trailing  edge  contour. 

It  was  found  that,  in  the  range  of  small  negative  angles  of  attack,  the 
drooped  trailing  edge  yields  better  aerodynamic  performance  than  the 
spiralled  one. 

2.  A  major  difficulty  in  using  the  viscous-inviscid  interaction 
method  along  with  the  optimization  code  is  its  inability  to  provide  fairly 
smooth  gradients  of  the  objective  function.  It  is  possible,  however,  to 
make  optimum  calculations  for  maximum  lift  by  reducing  the  number  of 
design  variables  to  the  minimum  and  removing  some  irregularities  in  the 
integration  of  aerodynamic  force  coefficients  from  the  surface  pressure 
distribution.  For  minimizing  the  drag,  further  improvement  of  the 
analysis  method  is  required. 

3.  Based  on  a  limited  comparison  of  the  theoretical  results  with 
experimental  data  obtained  with  negligible  wall  blockage  effect,  it  was 
found  that  the  Dvorak-Kind  method  overpredicts  the  lift  but  underpredicts 
the  drag. 
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Figure  4  -  Theoretical  and  Experimental  Pressure  Distributions 
Over  Model  103  Circulation  Control  Airfoil  at  ■  0.3  and 

a  »  -0.01  Degree 
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Figure  5  -  Program  Flow  Chart 


re  6  -  Overlays  of  the  Combined  Program 


Figure  8b  -  Lagrangian  Interpolations 
Figure  8  -  Integration  of  Pressure  Drag 
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Figure  11  -  Pressure  Drag  Coefficients  in  Two-Variable  Space 
at  =  0.3,  =  0.04,  and  a  =  -5  Degrees 
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SEARCH  POSITION 


Optimal"  Airfoil  Shape  Resulting  from  Maximizing  at  M£ 
C  =  0.04,  and  a  =  -2  Degrees 


TABLE  1  -  LIFT  AND  DRAG  COEFFICIENTS  OF  MODEL  103  CIRCULATION  CONTROL 
AIRFOIL  AT  =  0.3  AND  a  =  -0.01  DEGREE 


Run 

C. 

1-  • 

CIRCON  Results 

DTNSRDC  Experiment 
(7  ft  x  10  ft  Transonic 
Wind  Tunnel) 

n 

CD 

Up 

cDf 

CD 

cL 

CDp 

CD 

212 

0.0073 

0.6842 

0.0048 

0.0047 

0.0022 

0.4980 

0.0084 

0.0117 

213 

0.0105 

0.8641 

0.0079 

0.0048 

0.0021 

0.7008 

0.0148 

0.0208 

214 

0.0111 

0.8835 

0.0115 

0.0047 

0.0051 

0.7542 

0.0181 

0.0211 

215 

0.0160 

1.1695 

0.0067 

0.0043 

-0.0050 

1.0932 

0.0556 

0.0236 

216 

0.0220 

2.1620 

0.0227 

0.0044 

0.0050 

1.3634 

0.1112 

0.0413 

217 

0.0384 

3.1714 

0.0698 

0.0043 

0.0357 

1.9203 

0.1797 

0.1106 

218 

0.0291 

2.5010 

0.0308 

0.0046 

0.0063 

1.6575 

0.1714 

- 

NOTES:  1.  (CD)  =  CD  +  CD  -  C_ 

Theory  ^ 


2.  (Cn)  is  obtained  by  integrating  the  wake  rake  survev, 

u  Exp 

independent  of  surface  pressure  distribution. 

3.  C^  and  CD  versus  C;  are  plotted  in  Figures  2  and  3. 


TABLE  2  -  COORDINATES  OF  BASELINE  AIRFOILS 


Nijmmkr 

n 

M*FOILS» 

NAFAfL 

* 

3 

NUMHCS 

or 

X-STATIONS 

•  NX$TA 

m 

30 

X-L0CAT10NS  or  COORDINATES 

0.0000 

.00*0 

.0200 

.0400 

.0700 

•  1000 

•  1500 

.2000 

.2500 

.3250 

.*000 

•  5000 

,5750 

.6500 

.7000 

,7500 

.eooo 

.0500 

.8750 

.9000 

.9250 

,9500 

•  9600 

.9700 

.9000 

.9075 

.9925 

•  9965 

.9990 

1*0000 

BASIS  VICTORS 


VECTOR  NUMBER  1 

BASIS  AIRFOIL  A 

LOVER  SURFACE  COORDINATES 

.36O0SE-O2  -.IUB5E-01  -.23991F-01 
-.65097E-01  -.68802E-01  -.709T0E-01 
-.63903E-01  -.5961SE-01  -.S692TE-01 
-.30313E-01  -.2453SE-01  -.1917*5-01 

UPPER  SURFACE  COORDINATES 


■.33318E-CI 
•  •721 15E-r  1 
■.53708E-CI 
..U8S3E-CI 


-.4241SE-01 

-.71826E-01 

-.A953*E-01 

-.60935E-02 


■•48T9SE-01 

..70529E-01 

■•♦3A20E-01 

.36000E-02 


•.562TTE-01 

..6900SE-01 

•.*00*85-01 


•.61410E-0] 

•.6682SC-01 

-.351T2E-01 


.360B5E-02 

.11I88E.00 

.10828E-00 

.♦2280E-0I 

.itsabe-01 

, 1 228BE  *  0  0 
.9T481E-01 
.3A066E-01 

.33288E-0I 
• 1299SE.00 
.91070E-01 
•26963E-01 

.448I2E-01 

.13389E.C0 

.83862E-CI 

.19I80E-01 

•61805E-01 
• 132B8E.00 
•75162E-01 
•11501E-01 

.73719E-01 
• !2846E*00 
.43562E-01 
. 36000E-02 

.89544E-0I 

•12353E-00 

•57494E-01 

•  1 0197C-00 
. 1 1687E*00 
•50032C-01 

VECTOR  NUMBER  2 

BASIS  AIRFOIL  B 

LOWER  SUPFACE 
.36O09E-O2 
-.65097E-01 
-.63803E-01 
-.290B5E-01 

COORDINATES 
- , 1 1 485E-0 1  - 
-.68B02E-01  - 
-.59615E-01  - 
-.2E4T5E-01  - 

•2399IE-01 
•70978E-01 
•S6927E-01 
• 16590E-01 

-.3I318E-0I 

-.72115E-01 

-.5A0I9E-01 

•.9C999E-02 

-.♦2415E-01 
—  T1826E-01 
-.49806E-0I 
-.13898E-02 

-.48795E-01 
-.T0529E-01 
-.43395E-0 1 
•89000E-02 

-.542TTE-01 

-.49005E-01 

-.39933E-01 

-.61410C-01 
-.66828E-01 
—  35247E-01 

UPPER  SURFACE 
•360B5E-02 

•  me?E»oo 

• 1 O02flE*OO 
.49356E-01 

COORDINATES 

.IT5ABE-01 

.1228BE.00 

,9T»81E-01 

.A217TE-01 

.332B8E-01 
, 12995E. 00 
.91070E-01 
• 34969E-0 1 

,4fte.2E-0I 
,133  9E-00 
,83j6jE*01 
,2T5l*E-01 

•61805E-01 

.13288E.00 

.75692E-01 

. 19944E-01 

.73T19E-01 

.12B46E.00 

•66235E-01 

•89000E-02 

•  8954*f  — 11 
.12353E-1 ) 
•61809E-C1 

•  1 0197C*00 
•11607E-00 
.S63T3E-01 

VECTOR  NUMBER  3 

BASIS  AIRFOIL  C 

LOt*  TR  SURFACE 
•  *»6085E-02  . 

-  .A5097E-01 

-  6J0O3E-O1 

-  23928E-01 

COORDINATES 
-.U485E-01  - 
-.68802E-01  - 
-.S9615E-01  - 
-.16004E-01  - 

•23991E-01 

.70970E-01 

.569275-01 

.83T0SE-02 

-.333' EE-01 
-.721 i5c-0l 
-.53T87E-01 
..10980F.03 

-.42415C-01 
-.T1828E-01 
—.49401 E-0 1 
.80892E-Q2 

-.48T95E-01 
-.T0529C-01 
- . 4 1 828E-0 1 
.194771-01 

-.542775-01 

-.690052-01 

-.37491E-01 

-.61410E-01 

-.66828E-0I 

-•31649C-01 

U*PF*  SURFACE 
•30O05E-O2 

•  1 1 1 86E*00 

•  1 0828E*00 

COORDINATES 

.17548E-01 

.122S8E-00 

.9T481E-01 

.3328BE-01 
. I 2995E .00 
•9I070E-01 

.A8B12E-01 
. 1 J389E.00 
•83883E-01 

•6180SE-Q1 
• 13288E.00 
. 75669E-0 1 

.T3719C-01 

•12844S.00 

•44082E-01 

.89564E-01 
.123532. JO 
.617155*01 

.101971*00 
. 1148TE-00 
.S69S2E-01 

•51181E-01 

.45T43E-01 

.40807E-01 

, 3A6A7E-01 

. 27302E-01 

•1947TC-01 

THIS  PAGE  IS  BEST  <*UAim  PRACXlCAisiiA 
FROM  COPY  FUkSlSKEd.'  10  DDC 
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TABLE  3  -  RUN  MATRIX  FOR  CIRCULATION  CONTROL  AIRFOIL  OPTIMIZATION  BASED  ON  BASELINE 

SHAPES  Y, ,  Y. ,  AND  Y- 


-5  deg 


TABLE  3  (Continued) 
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